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Retardation of the Onset of Turbulence by Minor Viscosity Contrasts 
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Motivated by the large effect of turbulent drag reduction by minute concentrations of polymers 
we study the effects of minor viscosity contrasts on the stability of hydrodynamic flows. The key 
player is a localized region where the energy of fluctuations is produced by interactions with the 
mean flow (the "critical layer"). We show that a layer of weakly space-dependent viscosity placed 
near the critical layer can have very large stabilizing effect on hydrodynamic fluctuations, retarding 
significantly the onset of turbulence. The effect is not due to a modified dissipation (as is assumed in 
theories of drag reduction), but due to reduced energy intake from the mean flow to the fluctuations. 
We propose that similar physics act in turbulent drag reduction. 



The addition of small amounts of polymers to hydro- 
dynamic systems produces dramatic effects on phenom- 
ena such as the transition to turbulence, vortex forma- 
tion and turbulent transport [jlj. The phenomenon that 
attracted most attention was, for obvious reasons, the 
reduction of friction drag by up to 80% when a very 
small concentration of long-chain polymers were added 
to turbulent flows [^JD . In spite of the fact that the phe- 
nomenon is robust and the effect huge, there exists no ac- 
cepted theory that can claim quantitative agreement with 
the experimental facts. Moreover, it appears that there is 
no mechanistic explanation. In the current theory that is 
due to de Gennes |4|,|5| one expects the Kolmogorov cas- 
cade to be terminated at scales larger than Kolmogorov 
scale, leading somehow to an increased buffer layer thick- 
ness and reduced drag, but how this happens and what is 
the fate of the turbulent energy is not being made clear. 

In this Letter we propose that the crucial issue is in 
the production of energy of hydrodynamic fluctuations 
by their interaction with the mean flow. For the sake of 
concreteness we examine a simple laminar flow and its 
loss of stability, and show how small viscosity contrasts 
lead to an order of magnitude retardation in the onset 
of instability of "dangerous" disturbances. In this model 
everything is explicitly calculable, and we demonstrate 
that nothing special happens to the dissipation. Rather, 
it is the energy production that is dramatically reduced, 
giving rise to a large effect for a small cause. At the 
end of this Letter we argue that similar physics may very 
well be at the heart of turbulent drag reduction, but we 
stress that the phenomenon discussed below is interesting 
by itself and well warrants an experimental confirmation. 

It is well known that parallel Poiseuille flow loses its 
stability at some threshold Reynolds number Re=Reth 
(close to 5772). It is also well known that the instability 
is "convective" , with the most unstable mode having a 
phase velocity c. Analytically it has the form 

4>(x,y,t) = ±{(f>(y) exp[ik(x - ct)\ + c.c.} cxp( 7 i) , (1) 

where <p(x,y,t) is the disturbance streamfunction and 
4>{y) is the complex envelope of <j)(x, y, t) . We have chosen 
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FIG. 1. Schematic of the flow: the fluid near the walls has 
a viscosity /ii , and that flowing at the center is of viscosity [i? . 
In the mixed layer (of width q) the viscosity varies gradually 
between fii and /j.2- The parameter p controls the position of 
the mixed layer. 

x and fc as the streamwise coordinate and wavenumber of 
the disturbance, c as the phase speed and t as time. 7 is 
the increment of instability. What is not usually empha- 
sized is that the main interactions leading to the loss of 
stability occur in a sharply defined region in the channel, 
i.e. at a layer whose distance from the wall is such that 
the phase velocity c is comparable to the velocity of the 
mean flow. We refer to this layer as the "critical" layer. 
It is thus worthwhile to examine the effect on the stabil- 
ity of Poiseuille flow of a viscosity gradient placed in the 
vicinity of the critical layer. Following Q we examine 
a channel flow of two fluids with different viscosities Hi 
and H2, see Fig. ^ 

Observing that the inferred effective viscosity in poly- 
mer drag reduction increases towards the center by about 
30% over about a 1/3 of the half-channel 0, we choose 
fj,2 = 1 and m = = 0.9, with all the viscosity dif- 

ference of 0.1 concentrated in a "mixing" layer of width 
0.1, leading to comparable viscosity gradients. The ob- 
servation that we want to focus on is shown in Fig. 2: 
the threshold Reynolds number for the loss of stability 
of the mode as in Eq. (Q) depends crucially on the posi- 
tion of the mixing layer. When the latter hits the critical 
layer the threshold Reynolds number for the loss of sta- 
bility reaches as much as 88000. In other words, one can 
increase the threshold of instability for a given mode 15 
times, and by making the mixing layer thinner one can 
reach even higher threshold Reynolds values. 
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FIG. 2. The dependence of the threshold Reynolds number 
on the position of the viscosity stratified layer for m = 0.9. 
The dashed line pertains to the neat fluid. Note the huge 
increase in R t ^ within a small range. This occurs when the 
stratified layer overlaps the critical layer. 

In this Letter we analyze the physical origin of this 
huge sensitivity of the flow stability to the profile of the 
viscosity. The stability of this flow is governed by the 
modified Orr-Sommerfeld equation S 



ik [(</>* 



k 2 cf) (U-c-ij)- U"cf\ = J- 



^ (4) + 2 M y 



2k 2 fx'( l 



fcV) 



(2) 



in which <f>, U and p are functions of y. The boundary 
conditions are <^(±1) = 0'(±1) = 0. All quantities have 
been non-dimensionalised using the half-width H of the 
channel and the centerline velocity Uq as the length and 
velocity scales respectively. The Reynolds number is de- 
fined as Re = pUoH/^,2, where p is the density (equal 
for the two fluids). The primes stand for derivative with 
respect to y. At y = 0, we use the even symmetry condi- 
tions 4>(0) = 1 and 4>'(0) = 0, as the even mode is always 
more unstable than the odd. 

Since the flow is symmetric with respect to the channel 
centerline, we restrict our attention to the upper half- 
channel. Fluid 2 occupies the region < y < p. Fluid 1 
lies between p + q < y < 1. The region p < y < p + q con- 
tains mixed fluid. The viscosity is described by a steady 
function of y, scaled by the inner fluid viscosity p,2- 

H(y) = l, for 0<y<p, (3) 
[i(y) = l + (m~l)e[W-15t + 6e], 0<£<1, (4) 
fi(y) = m , for p + q < y < 1 , (5) 

(6) 

Here £ = (y—p) jq is the mixed layer coordinate. We have 
assumed a 5th-order polynomial profile for the viscosity 
in the mixed layer, whose coefficients maintain the 




FIG. 3. Profiles of the normalized viscosity p(y) and nor- 
malized velocity U(y) and the second derivative U (y) for 
m = 0.9 (solid lines) and m = 1.0 (dashed lines). The mixed 
layer is between the vertical dashed lines. 

viscosity and its first two derivatives continuous across 
the mixed layer. The exact form of the profile is unim- 
portant. For a plot of the profile m — 0.9, see Fig. [| 

The basic flow U(y) is obtained by requiring the ve- 
locity and all relevant derivatives to be continuous at the 
edges of the mixed layer: 

U(y) = l-Gy 2 /2, for y<p, (7) 
U(y) = U(p)-Gf^dy, for P <y< P + q, (8) 

U(y) = G (1 - y 2 ) /2m, for y > p + q . (9) 

Here G is the streamwise pressure gradient. 

It can be seen, comparing the mean profile U(y) to 
that of the neat fluid (cf. Fig. ||), that nothing dra- 
matic happens to this profile even when the mixing layer 
is chosen to overlap a typical critical layer. Accordingly 
we need to look for the origin of the large effect of Fig. 
2 in the energetics of the disturbances. To do so, re- 
call that the streamwise and normal components of the 
disturbance velocity u(x,y,t) and v(x,y,t) may be ex- 
pressed via streamfunction as usual: 

u(x,y,t) = d<p/dy and v{x, y, t) = — d<p/dx . (10) 

These functions may be written in terms of complex en- 
velopes similar to Eq. (|j): 



(x,y,t) = i{u(y)exp[ifc(iz:-c£)]+c.c.}exp(7i), (11) 
(x, y, t) = i{v(y) exp [ik(x — ct)] + c.c.} exp(7£) . 



The pressure disturbance p is defined similarly. 

Define now a disturbance of the density of the kinetic 
energy 



E(x, y,t) = | [u(x, y, t) 2 + v(x, y, tf 



(12) 



we can express the mean (over x) density of the kinetic 
energy as follows: 

E(y, t) = (E(x, y, f)^ = £ (y) exp (2 7 t) , (13) 

e(v) = l(Hv)\ 2 +*Wv)\ 2 ) ■ 

The physics of our phenomenon will be discussed in 
terms of the balance equation for the averaged distur- 
bance kinetic energy. Starting from the linearized Navier- 
Stokes equations for u and v, dotting it with the distur- 
bance velocity vector, averaging over one cycle in x and 
using Eqs. ©-(HI) leads to 

2 7 £ (y) = V • J(y) + W+ (y) - W„ (y) , (14) 

where the energy flux J(y) in the y direction, rates of 
energy production (by the mean flow) W+ (y) and energy 
dissipation (by the viscosity) W- (y) are given by 

T( n_ [u(y)p*(y) + c.c] l 

J ( y > = 4~ p + ^(yW£(y), (15) 

W+(y) = -\u'{y) [u(y)v*(y) + c.c] , (16) 
W^(y)^^-! [ 2k 2 £(y) + ^[\u'(y)\ 2 + kV(y)| 2 ]} . (17) 

The superscript * denotes the complex conjugate. To 
plot these functions we need to solve Eq. (||) as an eigen- 
value problem, to obtain c, 7, and (f>(y) at given Re and 
k. The value of c determines the position of the critical 
layer. 

It is convenient to compute and compare the space 
averaged production and dissipation terms T + and T_ 
defined by: 

r± = / W±(y)dy I f £(y)dy . (18) 
Jo 1 Jo 

The local production of energy can be positive or neg- 
ative, indicative of energy transfer from the mean flow 
to the disturbance and vice-versa respectively. The pro- 
duction in one region (where W+(y) > 0) can be partly 
canceled out by a "counter-production" in other region 
(where W+(y) < 0). 

The use of these measures can be exemplified with the 
neat fluid (m = 1.0 here). The laminar flow displays its 
first linear instability at a threshold Reynolds number 
of Reth — 5772, which means that the total production 
r + across the layer becomes equal to the total dissipa- 
tion r_ at this value of Re. Examining Fig. |] we can 
see that the disturbance kinetic energy is produced pre- 
dominantly within the critical layer, where the basic flow 
velocity is close to the phase speed of the disturbance, 
while most of the dissipation is in the wall layer. 
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FIG. 4. Energy balance: production W+(y), solid line; dis- 
sipation W-(y), dot-dashed line, Re = 5772. Top: m = 1, 
r+ = T_ = 0.0148. Bottom: m = 0.9, p = 0.3, F+ = 0.0158, 
r_ = 0.0148. In this and all the subsequent figures the solid 
vertical lines show the location of the critical lines, whereas 
the region between the dotted lines is the mixed layer. 



The balance is not changed significantly when the vis- 
cosity ratio is changed to 0.9 so long as the mixed layer 
is not close to the critical layer. There is a small region 
of production and one of counter-production within the 
mixed layer, whose effects cancel out, leaving the system 
close to marginal stability. 



We now turn our attention to Fig. g, in which the 
main point of this Letter is demonstrated. Here, the 
Reynolds number is the same as before, but the mixing 
layer has been moved close to the critical layer. It is im- 
mediately obvious that the earlier balance is destroyed. 
The counter-production peak in the mixed layer is much 
larger than before, making the flow more stable. The 
wavenumber used is that at which the flow is least stable 
for the given Reynolds number at this p. For m = 0.9, 
the threshold Reynolds number is 46400. Fig. ?? shows 
the energy balances at marginal stability - the picture is 
qualitatively the same here as at Re rs 5772 for the neat 
fluid. 
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FIG. 5. Energy balance for m = 0.9, p = 0.85. Energy re- 
duction W+(y), solid line; dissipation W-(y), dot-dashed line. 
Upper panel: stable flow at Re = 5772 (with T+ = -0.0114, 
r_ = 0.0122). Lower panel: Marginal flow at Re = 46400 
(with T+ = T_ = 0.0053). 



albeit at a much higher Reynolds number than the pri- 
mary mode; when we stabilize a given mode substantially, 
we should watch out for other pre-existing/newly desta- 
bilized modes which may now be the least stable. 

Finally, we connect our findings to the phenomenon of 
drag reduction in turbulent flows. Since the total dissi- 
pation can be computed just from the knowledge of the 
velocity profile at the walls, any amount of drag reduc- 
tion must be reflected by a corresponding reduction of 
the gradient at the walls. Concurrently, the energy in- 
take by the fluctuations from the mean flow should reduce 
as well. Indeed, the latter effect was measured in both 
experiments M and simulations pQ|,p!l| . The question is 
which is the chicken and which is the egg. In our calcula- 
tion we identified that the reduction in production comes 
first. From Figs. 4 and 5 (upper panel) which are at the 
same value of Re we see that the dissipation does not 
change at all when the mixing layer moves, but the pro- 
duction is strongly affected. Of course, at steady state 
the velocity gradient at the wall must adjust as shown 
in Fig. 5 (lower panel). We recognize that in a turbulent 
flow there are a number of modes that interact, but we 
propose that a similar mechanism operates for each mode 
at its critical layer, where both elastic and viscous effect 
determine the mean flow. In the present calculation we 
can consider all the putative unstable modes, and con- 
clude that with a viscosity gradient similar to that seen 
in polymeric turbulent flows the threshold Re goes up 
five times (to 31000). We leave the confirmation of this 
prediction to future experiments. 
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The main factor determining the stability is the pro- 
duction, which is driven by the phase change caused by 
the viscosity stratification. The dissipation on the other 
hand depends only on Reynolds number and does not 
respond disproportionately to changes in viscosity. In 
neat fluids, the term containing U"{y) in ^ is always 
of higher order within the critical layer. However, with 
the introduction of a viscosity gradient within the critical 
layer, the gradients of the basic velocity profile will scale 
according to the mixed layer coordinate £. An analysis 
in the critical layer indicates that for q < 0(i?e -1 / 3 ), 
the term containing U"/(U — c) is now among the most 
dominant. Any reasonable viscosity gradient of the right 
sign will pick up this term, leading to vastly enhanced 
stability. Note the dramatic effect in U" in Fig. 3. 

Indeed, in the light of this discussion we can expect 
that the large effect of retardation of the instability would 
even increase if we make the mixing layer thinner. This is 
indeed so. Nevertheless, one cannot conclude that insta- 
bility can be retarded at will, since other disturbances, 
differing from the primary mode, become unstable first, 
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